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NOTE  ON  SHOCK  REFLECTION  COEFFICIENT 

ABSTRACT 

The  shock  reflection  coefficient,  R^,  is  defined  here  as  the  ratio 
of  pressure  drop  across  a  reflected  disturbance  from  a  shock  wave  to  that 
across  the  incident  disturbance  wave  (which  is  generated  downstream  of 
the  shock).  The  treatment  of  shock  reflection  coefficients  in  the  work 
of  previous  investigators  is  reviewed  briefly.  A  derivation  is  given  of 
a  general  expression  for  R^  in  terms  of  shock  curvature  for  nonequili¬ 
brium  axisymmetric  flow.  Sample  computations  show  the  effects  of  non¬ 
equilibrium  flow  and  three-dimensionality  on  the  behavior  of  R  .  It  is 

,  w 

shown  that  R^  can  become  infinite  in  these  situations,  although  for  two- 
dimensional  ideal  gas  flow  its  magnitude  is  always  finite  and  usually 
very  small. 
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LIST  OF  SYMBOLS 

f  Nonequilibrium  production  term,  Eg..  (4) 

K  Shock  wave  curvature 
w 

M  Local  Mach  number 

n  Arc  length  on  curve  normal  to  streamlines 
p  Pressure 
q,  Flow  speed 

Shock  reflection  coefficient 
s  Arc  length  on  streamline 
T  Temperature 

» 

x  Coordinate  in  direction  of  free  stream 
y  Coordinate  normal  to  free  stream 
P  Shock  wave  angle 
y  Ratio  of  specific  heats 

e  -  0  for  two-dimensional  flow;  =  1  for  axisymmetric  flow 
7)  Are  length  on  left  running  Mach  line 
0  Direction  of  velocity  vector 
X  =  p-0 

p  Mach  angle,  =  arc  sin  (l/M) 

§  Arc  length  on  right  running  Mach  line 
p  Density 
Subscripts : 

®  Free  stream  value 
w  Value  at  shock  wave 


7 


INTRODUCTION 


The  shock  reflection  coefficient,  designated  here  by  Rw,  is  a 
quantity  of  interest  in  the  study  of  interaction  of  shock  and  Mach  waves, 
particularly  in  connection  with  application  of  the  shock-expansion  method 
to  computing  supersonic  flows.  Essentially  it  provides  a  measure  of  the 
strength  of  a  reflected  disturbance  from  a  shock  wave  relative  to  the 
strength  of  the  incident  disturbance  wave  (which  is  generated  downstream 
of  the  shock). 

12  3* 

Lighthill  ,  Chu  ,  and  Chernyi  have  obtained  the  reflection 
coefficient  for  two-dimensional  ideal  gas  flow.  (See  page  180  of 
Reference  5  for  a  brief  discussion  of  errors  in  References  1  and  2.) 
Graphs  of  computations  found  in  References  5  and  k  show  that  the  magni¬ 
tude  of  is  small  except  near  the  region  of  shock  detachment,  and 
that  |rJ<1  for  finite  Mach  number.  Note,  however,  that  R^  =  -  1  for  a 

straight  shock  and  finite  flow  deflection,  according  to  the  definition 
in  Equation  (l).  Eggers  and  Syvertson^  consider  a  closely  related 
quantity  called  "disturbance  strength  ratio.”  Their  calculation  for  a 
case  of  air  in  vibrational  equilibrium  (l^  =  10,  T^  =  278°  K)  indicates 
that  caloric  imperfections  in  the  gas  still  leave  the  disturbance 
strength  ratio  small  compared  to  unity  if  y  does  not  decrease  appreciably 
below  1.5* 

We  now  wish  to  investigate  nonequilibrium  effects  on  the  shock 
reflection  coefficient;  here  the  shock  curvature  appears  explicitly.  In 
addition,  we  wish  to  examine  the  effect  of  three-dimensionality,  after 
extending  the  definition  of  R^  to  axisymmetric  flow. 

EQUATIONS 

Refer  to  Figure  1;  the  shock  reflection  coefficient  is  defined  as 
follows  (as  on  page  252  of  Reference  6): 

lin  pc  •  %  .  lia  (ap/an)B  ^ 

IV  -  A5,4H>  pb  -  pa  ■  idp/d5)A  ^  > 

*  Superscript  numbers  denote  references  which  may  be  found  on  page  15. 
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OBLIQUE  SHOCK 


(Xa/9-0).  ACABb(m+X),  Z-ACB«(/t-X) 
A~B«A£,  BC» Aij 

- »-x 


FIGURE  1.  CROSS-SECTIONAL  VIEW  OF  FLOW  FIELD 


Here  the 
positive 
get 


By  means 

dp/dTl  =  (dp/ds)  cos  p,  +  (dp/dn)  sin  p 

dp/dj-  =  (dp/ds)  cos  p  -  (dp/dn)  sin  p 
and  the  following  relation  valid  at  the  shock  wave: 

(d/dn)w  =  (esc  x)w  [K^  d/dg  -  (cos  X)  d/ds]w 
the  shock  reflection  coefficient  is  expressed  as 

V  Sln<*-l>[(sin,:§K’  -  isin(^)J 
and  dp/ds  is  given  ty  * 

D  to  =  (KV  cso  C  ip  cot  x  +  PI2  If-)  +  *  PI2  sin  0  +  pqSf  (4) 

where  D  =  1  -  M2  +  cot2!,  and  f  is  the  nonequilibrium  ’’production  term" 

(f  =  0  for  frozen  or  equilibrium  flow);  e=  0  or  1  for  two-dimensional  or 
axisymmetric  flow,  respectively. 

CALCULATIONS 

We  consider  dissociating  air,  using  the  model  presented  in  Reference 
7,  as  an  example  of  nonequilibrium  effects  on  R^j  for  this  flow  the  air 
is  in  vibrational  equilibrium  at  the  shock.  A  typical  set  of  results  is 
sho»m  in  Figure  2  for  a  wide  range  of  shock  curvature.  It  is  seen  that 
for  each  value  of  there  is  a  critical  deflection  angle  producing 
infinite  R^;  this  results  from  the  vanishing  of  dp/d§  in  Equation  (2). 


positive  5  and  T]  directions -make  angles  of  p(<  90°)  with  the 
streamline  direction.  Applying  the  law  of  sines  to  A  ABC,  we 

of  the  directional  derivative  relations 


This  formula  ie  found  in  section  (5  of  Reference  7t  where  a  typographical 
error  occurs ;  the  multiplication  sign  should  he  replaced  by  a  plus  sign. 
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SHOCK  REFLECTION  COEFFICIENT  FOR  NONEQUILIBRIUM  DISSOCIATING 
AIR,  TWO  DIMENSIONAL  FLOW 


For  deflections  smaller  than  this  critical  value  the  behavior  of  Rw  is 
similar  to  that  for  ideal  gas  flow  except  for  the  sharp  rise  at  the 
asymptote;  beyond  the  asymptote  R^.  approaches  -  1  at  the  detachment  angle. 

When  e  is  taken  equal  to  1  in  Equation  (4),  the  definition  of  shock 
reflection  coefficient  [ (Equation  (l)]  can  be  applied  to  axisymmetric 
flow.  Figure  3  shows  some  typical  results  in  ideal  gas  for  several 
values  of  the  parameter  y!^.  (Generally  for  flows  over  smooth  ogival 
bodies  of  revolution,  -0.1  <  *0.)  The  behavior  of  R^  here  differs 

considerably  from  that  for  two-dimensional  flow.  There  is  a  critical 
deflection  angle  for  each  yK^  and  a  limiting  Rv  value  of  -  1  as 
detachment  angle  is  approached.  Formally,  yK^.  ■  -  00  corresponds  to  two- 
dimensional  flow,  and  the  figure  shows  the  approach  of  the  curves  with 
decreasing  yK^  to  the  two-dimensional  one. 
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